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The introduction of hquid '^He into siUca aerogel provides us a with model system in which to 
study the effects of disorder on the properties of a strongly correlated Fermi liquid. The transport 
of heat, mass and spin exhibits cross-over behavior from a high temperature regime, where inelastic 
scattering dominates, to a low temperature regime dominated by elastic scattering off the aerogel. 
We report exact and approximate solutions to the Boltzmann-Landau transport equation for the 
thermal conductivity of liquid '^He, including elastic scattering of quasiparticles by the aerogel and 
inelastic quasiparticle collisions. These results provide quantitative predictions for the transport 
properties of liquid ''He in aerogel over a wide range of pressure, temperature and aerogel density. 
In particular, we obtain a scaling function, _F(T/T*), for the normalized thermal conductivity, k/kci, 
in terms of a reduced temperature, T/T*, where T* is a cross-over temperature defined by the elastic 
and inelastic collision rates. Theoretical results are compared with the available experimental data 
for the thermal conductivity. 
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I. INTRODUCTION 



■ Aerogels are extremely low density solids formed as a rigid 

■ network of silica strands and clusters having a typical diame- 
iter of 30-50 A and porosities [g) above 99% They turn out 
' to be an ideal system for studying the effects of quenched ran- 
,dom disorder on the otherwise pure, ordered phases of liquid 
'■^He. When impregnated with liquid •^Hc, the aerogel is found 
\ to have dramatic effects on the transport properties and the 

■ phase diagram of liquid ^He, although basic thermodynamic 
'features characteristic of the Fermi liquid, such as the com- 
. pressibility, magnetization and heat capacity, are essentially 

■ unchanged The low temperature transport of mass, heat 
,and magnetization is substantially reducedf^i^i^ In addi- 
■tion, the superfluid transition temperature as well as the 
Isuperfluid order parameter are strongly suppressed relative 

■ to their bulk values i^i^ 

In this paper we consider the effects of scattering of ^He 

■ quasiparticles off a uniformly distributed random potential 
representing the aerogel structure, referred to as the "homo- 
geneous scattering model" (IISM)i^ We obtain exact and 
approximate solutions to the Boltzmann-Landau transport 
equation for the thermal conductivity of liquid "^He, includ- 
ing both inelastic collisions between quasiparticles and elas- 
tic scattering of quasiparticles by the random potential. The 
chief inadequacy of the HSM is its neglect of the inhomoge- 
neous void-structure of the aerogel, or more generally meso- 
scopic correlations that are observed in static structure fac- 
tor, and to which the superfluid transition temperature is 
sensitivei22i2ii^ However, the transport properties are lim- 
ited by the mean free path for quasiparticles propagating 



ballistically within the aerogel, and hence are expected to be 
well accounted for in the framework of the HSM since the 
geometric mfp (£) is typically much longer than the aero- 
gel correlation length {^a), e.g. ^/^^ ~ 3 for a 98% aerogel. 
Possible corrections to transport processes resulting from a 
small distribution of large voids or to fractal correlations on 
mesoscopic length scales < within the aerogel arc not in- 
cluded in the analysis presented in this work. However, the 
exact solution to the transport equation for the two-channel 
scattering model discussed in this paper should be of value in 
identifying observable corrections associated with correlated 
disorder or corrections to the two channel scattering theory. 

Liquid '^He is a dense quantum liquid in which the interac- 
tions between Fermionic excitations (quasiparticles) are one 
to two orders of magnitude larger than the mean kinetic en- 
ergy per particle. These interactions lead to strongly renor- 
malized branches of Fermionic excitations, reflecting the cor- 
related motion of many '^He atoms, and the emergence of 
Bosonic excitations. The Fermionic excitations bear resem- 
blance to "^Hc atoms only in terms of their quantum numbers 
for spin (s = ±1/2) and fermion number (e ~ ±1)- The 
Bosonic excitations come with and without spin and can be 
understood in terms of pairs of Fermionic excitations, e.g. 
the phonons of zero sound. Finally, the coupling between the 
Bosonic and Fermionic excitations leads to finite lifetimes for 
both types of excitations<^i^"— 

Interactions between "^He quasiparticles enhance the collision 
rate for quasiparticles near the Fermi surface, leading to a 
significant reduction in the lifetime of a quasiparticle at the 
Fermi surface. However, Fermi statistics rescues the low- 
energy quasiparticles (as well as the Bosonic modes) 4S At 
low temperatures, T <^ Ef/kB « 1 K, the number of exci- 
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tations is low, n^p « {JtBT / Ef)n. Similarly, binary collision 
processes are confined to a small region of phase space near 
the Fermi surface, Ap « {UbT / Ef)pf . As a result the Pauli 
exclusion effect suppresses the quasiparticle collision rate?^ 



1 

Tin 



{W) {ksTf 



(1) 



where m* is the effective mass of a quasiparticle and (W) is 
the square of the transition matrix element for binary colli- 
sions averaged over the Fermi surface. 

Furthermore, since the density of excitations is low, n^p <C n, 
transport coefficients are given by formulae familiar from gas 
kinetic theory. In particular, the transport of heat is dom- 
inated by thermally excited quasiparticles with the thermal 
conductivity given by^ 



1 



K= ^i.CvVf){vfT^) , 



(2) 




P [bar] 



where c,. 



^^Nfkj^T is the low-temperature specific heat, 
Nf is the quasiparticle density of states at the Fermi energy, 
Vf is the Fermi velocity and VfT,^ is the transport mean- 
free-path for heat conduction, with ^ Ti„ ^ T~^. Thus, 
heat transport becomes very efficient in pure "^He with k ~ 

In aerogel, elastic collisions of ■^He quasiparticles with the sil- 
ica strands lead to a temperature-independent contribution 
to the mean free path and hence the quasiparticle scatter- 
ing rate. Thus, at sufficiently low temperatures, transport 
currents are limited by elastic scattering from the aerogel, 
whereas inelastic scattering of quasiparticles dominates at 
high temperatures. There is an intermediate regime where 
both mechanisms are important. The cross-over tempera- 
ture separating these regimes is estimated from the inelastic 
collision rate in the pure "^He in Eq. ([T]) and the mfp of the 
aerogel, £, which provides an estimate for the clastic collision 
rate, 
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constant . 



(3) 



Estimating from Eqs. ([T]) and ([3]) gives. 



FIG. 1: Cross-over temperature vs. pressure for aerogels charac- 
terized by the mfp as defined by Eqs. 11151 The curves correspond 
to ^ = 100 — 250 nm spaced by 10 nm. The red curve is for £ = 170 
nm (g ~ 98%). The cross-over temperature is obtained from Eq. 
[S]with (W) given by Eqs. (|174m77|) . Table [T] (see the appendix 
in Sec. Wl, and the relevant Fermi-liquid parameters from Refs. 

Eaiii. 



The cross-over from the high-tempcraturc regime dominated 
by inelastic scattering to the low-temperature regime domi- 
nated by elastic scattering off the aerogel is characteristic of 
all transport processes for ^He in aerogeh ^^i^^'^^i'^^ Below we 
report theoretical results for the heat transport coefficient of 
liquid '^He-aerogel, in the normal state, over the broad tem- 
perature range, Tc < T <^ Ef/kB- Since the bulk properties 
of "^He are well known, measurements of the transport coef- 
ficients can provide quantitative information on the effects 
of disorder on the transport of "^He quasiparticles through 
aerogel. 



II. TRANSPORT THEORY 



T <T, 
T>T,. 



(4) 



The cross-over temperature, T*, defined by ri„(T*) = Toi, is 
given by2^ 



8Ef/kB 
^/irkfi {W) 



(5) 



where (W) is the dimensionless quasiparticle transition prob- 
ability averaged over the Fermi surface (Eq. 11711 of the ap- 
pendix). For 98% porosity aerogel we estimate £ w 1700 A)2^, 
and thus T^, « 18 mK at p = 15 bar (see Fig. [1]). 



In normal liquid '^He at low temperatures the transport of 
mass, energy and magnetization is carried predominantly 
by fermionic quasiparticles, whose distribution in phase 
space, ripo-, is governed by the Boltzmann-Landau transport 
equation;^ 



dup 
dt 



Ip. (6) 



where p ~ (p, cr) denotes the momentum and spin of the 
quasiparticles. For spin-independent transport the quasipar- 
ticle energy 



ep = Ep + <5£p(r,i) + u<,,t(p,r,t) , 



(7) 
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is the sum of the equilibrium excitation energy, ep, the cou- 
phng to an external scalar or vector potential, M„^t(p, r, i), 
and the Landau molecular field energy. The latter arises 
from the interaction of a quasiparticle with the distribution 
of non- equilibrium quasiparticles, 

p'a' 

where 6np is the deviation of the distribution function from 
(global) equilibrium, 

Sup = rip - no(ep) . (9) 

For small disturbances from equilibrium the derivative of the 
equilibrium distribution function, 

confines the excitations to states that lie near the Fermi sur- 
face, epj, = fi. The interaction energy between two quasipar- 
ticles is given by /p,p', and in contrast to the distribution 
function, varies slowly with |p| in the vicinity of the Fermi 
surface. We can typically evaluate /p,p', as well as the den- 
sity of states, iV(e), on the Fermi surface, i.e. N{e) ~ Nj, 
and 2Nff{pfp,pfp') = ^^(p • p') = ^> ^Hp ' P')- The 
latter equality defines the dimensionless Landau parameters. 



A. Collision Integrals 

The right side of the transport equation, Jp, represents the 
change in the distribution function resulting from collision 
processes. We consider two scattering processes for '^He in 
aerogel: (i) elastic collisions of quasiparticles with "impuri- 
ties" representing the aerogel strands, and (ii) inelastic colli- 
sions between quasiparticles. The development of the trans- 
port theory for '^He-aerogel presented below, particularly the 
reduction of the transport equation in the low temperature 
limit, parallels that development by Baym and Pethick in 
their review on transport in pure liquid ^Hc;^ and extends 
Brooker and Sykes' work on the transport coefficients of 
Fermi liquids 

In our case the effects of aerogel scattering enter through a 
contribution to the collision integral. For quasiparticle scat- 
tering by the aerogel strands 

2 

Ip[ = -y'^w{pi,p2)S{ei -e2)Ki -np^], (11) 

where 'w(j)i,p2) is the transition rate for scattering of quasi- 
particles by the aerogel. 

For inelastic quasiparticle-quasiparticle collisions at low tem- 
peratures, ksT <^ Ef, only binary collisions are important. 
We denote t{pi,p2',P3,P4) as the scattering amplitude for 



binary collisions between quasiparticles with momenta and 
spin Pi = (pi, (Tj ). The labels pi and p2 refer to initial states 
while p3 and p4 refer to final states. Fermi's Golden Rule for 
the transition rate {pi,P2) ~^ (PSTPi) is: 

r = ^ \t{pi,p2;p3,P4)\'^ S{ei + £2 - £3 - £4) ■ (12) 

For a translationally invariant system with spin-rotation 
invariant interactions between quasiparticles the transition 
rate includes momentum- and spin-conserving delta func- 
tions, 

r = W{pi,P2;P3,P4)Sp 1-|-P2,P3-|-P4 '^fTl+(T2,fT3-|-(T4 

X (5(ei -f £2 - £3 - £4) ■ (13) 
where is a smooth function of p^ . 

The collision integral for binary scattering includes the phase 
space factors for collisions that both increase and decrease 
the population of the state pi (scattering " in" and scattering 
"out"). In particular, 

4"=- r(pi,p2;P3,P4) X (14) 

P2,P3,P4 

Kl«P2(l - "P3)(l - "P4) - (1 - "Pl)(l - "P2)"P3"-P4] ■ 

The sum over final states (p3,P4) is restricted to avoid double 
counting of equivalent states of identical particles related by 
exchange of 0^4. 

The collision integral vanishes when evaluated with a local 
equilibrium distribution function, i.e. /[rip*^ ] = 0. For the 
clastic scattering contribution to the collision integral (Eq. 
[TT|) this identity is obvious. For the inelastic collision integral 
it is less so, but follows from the identity, 

'5(£pi +£p2 -£p3 -£p4) X (15) 

'^o(£p3)"o(£p4)(l - '^o(£pi))(l - '^o(£p2))] = 0, 

where no{e) = l/(e^^'^^'^''-|-l) is the Fermi distribution. This 
identity is a consequence of local equilibrium and the con- 
dition of detailed balance between the scattering "in" and 
scattering "out" contributions to the collision rate. 

Although translational invariance is violated by the pres- 
ence of the aerogel medium, the aerogel is sufficiently dilute 
that the scattering rate by the aerogel impurities is typically 
small compared to excitation energies in the normal state, i.e. 
ft/r^ ^ ksTi^ In this limit the effects of aerogel scattering 
on the intermediate states that enter the inelastic collision 
integral can be neglected. Thus, momentum conservation 
holds for the binary collision integral for normal '^Hc in high 
porosity aerogels. At lower temperatures, e.g. in the super- 
fluid phase, or for lower porosity aerogels, this approximation 
breaks down. This limit requires a microscopic treatment 
of the effects of aerogel scattering on inelastic collision pro- 
cesses which is outside the scope of the phcnomcnological 
Boltzmann-Landau transport theory. 
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B. Linearized Transport Equation 



The transport coefficients of liquid '^He in arcogel are cal- 
culated from solutions of the linearized transport equation 
in steady-state. The particular solution depends on the 
nonequilibrium conditions that are established. For small 
deviations from equilibrium the nonequilibrium steady state 
is specified by a local equilibrium distribution function, 



,/.e. 



1 



1 +e(ep(r)-M(r))/feBT(r) ' 



(16) 



parametrized by a local temperature, T(r), chemical poten- 
tial, fj,(r) and quasiparticle energy, ep(r). The transport 
equation naturally separates by expanding about this local 
equilibrium distribution. 



6n^ 



J.e. 



(17) 



because the collision integral vanishes under the conditions 
of local equilibrium, /p[?ip'^ ] = 0. The linearized transport 
equation becomes, 



J.e. 



dn: 



I.e. 



Vp • Vr£p = SIp[Snp] 



(18) 



Similarly, the quasiparticle heat current is given by the trans- 
port of excitations with energy, $p = Ep — 



(21) 



D. Elastic Scattering Limit 



Transport properties of normal '^He in aerogel at sufficiently 
low temperatures, i.e. T <C T*, are limited by elastic scatter- 
ing of quasiparticles by the aerogel structure. In this limit 
the transport equation is given by Eq. ITOl) with the collision 



term of Eq. pT|) . The integral equation for Sup is 



/ dno 



Vp.Z(ep) = -5] z«(p,p')'5(ep-ep') 

p'cr' 

X [Sfip^ ~ (5np'^/] , (22) 



where 



Z(Cp) = 



VT + V^t 



(23) 



The left side of Eq. ([181) supplies the driving terms, e.g. V/i 
and VT, for the collision terms on the right side that act to 
restore equilibrium. Using Eq. (|16p the linearized transport 
equation reduces to 



9no 



T 



VT + V^t 



= «p[<5np], (19) 



where (5/p is the collision integral to linear order in bn^. 



The solutions to Eq. ([2^ are determined by the energy 
and momentum dependences of the driving term and are 
familiar from the theory of electron-impurity scattering in 
metals.— To proceed further we need the '^He quasiparticle- 
aerogel scattering probability, w{p,p'). 



E. Elastic Scattering Model 



C. Quasiparticle Currents 

The mass and heat currents are determined by the solution 
for Sfip of Eq. (|19p . In particular, the mass current is given 

by 

jm = ^ m*Vp Sfipa , (20) 

per 

where m* ~ pf/vf is the quasiparticle effective mass. This 
form for the mass current is applicable to interacting Fermi 
liquids which are Galilean invariant)^ In pure liquid ^He 
quasiparticle-quasiparticle interactions which give rise to the 
enhancement of the Fermionic mass are Galilean invariant. 
For liquid '^He-aerogel Galilean invariance is violated by 
quasiparticle scattering off the aerogel. However, the non- 
Galilean contribution to the effective mass is of the order of 
concentration of scattering centers, Us/n <C 1, and thus neg- 
ligible compared to the quasiparticle-quasiparticle effective 
mass enhancement. 



We model the aerogel as a distribution of local scattering 
centers represented by the potential, U{r) — "(i* ~ I^j)- 
The terms u(r — Hi) represent the potential provided by 
the aerogel scattering centers at the fixed positions, {Ri|i = 
l...A'^s}- For a random distribution of uncorrelated scattering 
centers the rate is proportional to the mean number density 
of scattering centers, = Ng/V. In the Born approximation 
the transition rate is related to the matrix elements of u, 

wipa;p'a') ^ n,^\{p'a' \u\pa)\\ (24) 

For stronger scattering the potential u is replaced by the 
i-matrix for quasiparticle scattering by aerogel strandsi^ 

We shall assume that the scattering by the aerogel is non- 
magnetic. This should be sufficient for describing transport 
processes in zero field, particularly if the aerogel strands are 
"coated" with a layer of solid *He. However, it is known 
that ^He atoms form a highly polarizable solid layer on the 
surface of the aerogel strands and that these nuclear spins ex- 
hibit a Curie-like spin susceptibility)^ Thus, spin-exchange 
scattering of ^Hc quasiparticles by localized and polarizable 
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^He spins may be relevant to magnetic transport processes 
and transport in relatively low magnetic fields. The simplest 
scattering model for ^He-aerogel assumes the ^He quasipar- 
ticles interact with the aerogel via an isotropic scattering 
potential)^ There is no preferred direction within the aero- 
gel and the scattering probability depends on the relative 
orientation of the initial and final quasiparticle momenta. In 
this case, 



Li(p,p') = w{p ■ p';^p) = ^wi{^p)ri{p ■ p') 



(25) 



where w;(Cp) is the scattering probability for quasiparticles 
with relative orbital angular momentum I, and Vi{x) is the 
corresponding Legendre polynomial. Note that |p| = |p'| for 
elastic scattering, and we have parametrized the functional 
dependence of w(p, p') on |p| by the energy = w/(|p| —Pf) 
measured relative to the Fermi surface. The probabilities for 
scattering in the orbital channels are proportional to, 



+1 



dx 



(26) 



(2; + i)y_i 2 

which vary smoothly with on the scale oi Ef. 

For an isotropic scattering medium we make the ansatz, 

5np = (^Ig) [vpr„,(Cp)] • Z(^p) . (27) 

The momentum sum is represented as 



p' 



(...) 



An 



(28) 



The terms (Srio/^Cp) and S{^p — £,p') confine |p'| — |p| 
Pf + ^p/vf, so we obtain 



= 2iV(ep)(«;(p-p';ep)(l-p-p')>P' 



(29) 



= / ~^(- • ■)• We introduce the scattering rate 



where (. . 

for orbital channel / 

1 _ 27V(Cp) dx 



fj^ -w{x;^p)Vi{x)^2N{^p)wi{^p) 



rii^p) (2Z 

(30) 

and express the transport scattering rate in terms of the I ~ 
and 1 = 1 (s- and p-wavc) scattering rates. 



1 



1 



1 



r^Cp) ro(ep) MCp) 



(31) 



Note that the density of states, N{^p), and the scattering 
probabilities, w/(Cp), vary slowly with excitation energy on 
the scale of the Fermi energy. Thus, for many cases of interest 
we can safely neglect the energy dependence of the scattering 
rate and evaluate T„i(5p) — t„i(0) = r„i. 



In particular, the mass current induced by a pressure gradi- 
ent at constant temperature is determined by the quasipar- 
ticle mobility defined by j,-„ = —i/'Vfi. For fcsT <C Ef Eq. 
(|10p is sharply peaked at the Fermi energy, and the quasi- 
particle mobility calculated from Eqs. (pO|) and p7| is to 
leading order in T/Ef, 



V= o^/P/("/^ol), 



(32) 



which has the intuitive interpretation of transport of momen- 
tum Pf over a distance of order the mfp, i = VfT^i, within 



the aerogel. Similarly, the heat current, Eq. (PT|) . induced bv 
an temperature gradient, j, = —k^T, defines the thermal 
conductivity, which to leading order in T/Ef is given by. 



27r2 



kBNf{vfkBT)e, 



(33) 



and also has the simple interpretation as the flux of thermal 
energy VfksT transported over a distance of order I. 



F. Two Channel Collision Integral 



For higher temperatures, i.e. T > T^, both elastic and inelas- 
tic collision processes limit transport currents. The transport 
coefficients are then calculated from 



dng 



T 



VT + V/i 



(34) 



where the right side of Eq. (|34p contains the linearized col- 
lision integrals for both inelastic and elastic scattering. The 
elastic collision integral follows immediately from Eq. ([TT|) , 



?pj- (35) 



Since the driving terms in Eq. p4p arc confined to excitation 
energies within fc^T of the Fermi energy we express 



Srir, 



dnp 



i^p. 



(36) 



where ipp^ measures the deviation of the equlibrium distribu- 
tion for quasiparticles with excitation energy ^p. at the point 
Pi on the Fermi surface. 

For l^pj^l ^ Ef energy conservation, combined with the 
phase-space restriction required by the Fermi distribution 
factors in Eq. ([55]). forces the scattered excitation energy to 
be confined to the low-energy shell, i.e. l^pj < ksT. Thus, 
slowly varying functions of Pi can be be evaluated with mo- 
menta, |pi| = IP2I ~ Pf + £,pi/vf, i.e. in close vicinity of 
the Fermi momentum. The scattering rate w reduces to a 
function of the directions of the momenta for incident and 
scattered excitations with excitation energies near the Fermi 
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energy, w(pi,p2) w;(Pi, P2; CpJ- For the isotropic scat- 
tering model and a driving term proportional to Vp^ • VT we 
set 

^P. = Vp. • VT ^(^pj , (37) 
in which case the elastic collision integral becomes, 



61- = ^^Hi^ no(epJ(l - no(CpJ) ^ 



T 



ol(CpJ 



^(CpJ. (38) 



Linearizing the inelastic collision rate in Eq. (|14p , and mak- 
ing use of Eq. yields. 



P2,P3,P4, 

X [M(pi)M^P2)i^ - '^o(Cp3))(i - Mipi))] 

[tppi -1pp3 -■^P4] ■ 



(39) 



Here we consider an un-polarized Fermi liquid in which the 
only spin-dependent interactions are those that arise from 
exchange symmetry. In this case, S,p. = S^p- and the distribu- 
tion functions arc independent of (7^. Thus, we can carry out 
the sum over the spin states (72.34. We can also eliminate 
one of the momentum sums, resulting in, 

2 

^^Pl W^(Pl,P2;P3,P4)'5(^pi +Cp2 -Cp3 -Cp4) 

P2.P3 

X K(?pi)n-o(^p2)(l - '^o(?p3))(l - "o(Cp4))] 

X [■0P1 + V'P2 - V'P3 - V'pi] > 



where P4 = Pi -I- P2 — Pa and 

1 



W 



n 



(40) 



(41) 



is the spin-averaged scattering rate; W^^- is the scattering 
rate for cti = (T2 =t and W^^ is the rate for cti = — 0-2 =t- 
The weight factors take into account the restriction to avoid 
double counting of equivalent states, so the remaining mo- 
mentum sums over p3 and P4 are unrestricted.— 



For l^pj I ^ Ef the energy and momentum conservation laws, 
combined with the phase-space restrictions required by the 
Fermi distribution factors in Eq. (|40|) . force all excitation 
energies to be confined to the low-energy shell, i.e. |^pj < 
ksT. In this limit slowly varying functions of Pi can be 
evaluated on the Fermi surface. Thus, the scattering rate 
W becomes a function of the directions of the momenta for 
quasiparticles on the Fermi surface, 



W^(Pl,P2;P3,P4) W^(Pl,P2;P3,P4) , 

and the momentum sums can be approximated by 



p 



(...)- Nj 



dfip, 
47r 



d^p, (...)• 



(42) 



(43) 



To carry out the angular integrations we adopt Abrikosov 
and Khalatnikov's paramctrization^ of W in terms of the 
angle 9 between the two incoming momenta, and (/), the scat- 
tering angle between the planes defined by n = pi x p2 and 
n' = p3 X p4. The integration over the direction ps is ex- 
pressed in terms of angles relative to the conserved direction 
of the total momentum, P = pi -I- p2, 



drip 

47r 



-(...) 



2tt 



+1 



d{cos a) 



(...), 



(44) 



Since p4 — y/ P'^ + p'^ — 2P ■ we have dp4 — 
{ps/pijP d{cos a), and for momenta near the Fermi surface, 
P ~ 2pfCos{e/2). Thus, 



d(cosa) 



P4 



2vfPf cos{e/2) 



(45) 



Also, the azimuthal angle (pj, is the scattering angle up to a 
fixed but arbitrary constant, thus, d4>3 = d(p. In the case of 
the integration over the incoming momentum direction p2 we 
choose the remaining momentum direction pi as the polar 
axis, 



d<j>2 f^^ d{cose) 



27r 



(...) 



(46) 



The binary collision integral then reduces to 



T f\2vfpf 



j c?Cp2 j c?Cp3 j d£.pj{(.p, +Cp2 - Cp3 - ^pJ ["■o(^pi)'^o(^p2)(l - »^o(^p3))(l - noi^pM 



X 



r+i d(cos0) /-^^ d0 f W{e,c^) \ [^^ dcp2,, ^ , , , , 



(47) 



I 

The inhomogcneous terms of the linearized transport equa- excitation energy, ^p^ — > — Cpi > a-ii^d momentum direction, 
tion dictate the symmetry of the solution with respect to pi. We separate the angular and energy dependences of 
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the non-equilibrium distribution function with the ansatz, 



(48) 
(49) 



where ^''^''(Cpi) ~ ^'^'"^H^^pi)- We can now carry out the 
integration over 02 for each term in Eq. (|47l) . 



2n 



Vfiip^+H^pJ + VT^'-'i^P,) , (50) 



where the direction cosines, Xi = pi ■ pi for i = 2,3,4, are 
simply related to {0,4'). The angular integrations decouple 
from the energy integrations which are confined to the low- 
energy shell near the Fermi surface. We define the average 
scattering rate 



{W) 



+1 



d(cos 61) r'' / w{e, 



27r V2cos(6'/2) 



(51) 



as well as the weighted averages, r.i = {xiW) / {W) . Changing 
^P3,P4 ~^ ~Cp3,P4 resulting collision integral reduces to 



1 N] 



T pfVf 



(W^>vp, • X no(^pJ X /(epj^(^'(epj +^^^^/ d^P.M^P.)K{Cp,+Cp,)ip^^H^P. 



, (52) 



where 



K{0 = 



no 



d^3 j dUSiC + (s+U)noi^3)no{U) 



d6rio(6)A'(C + 6) 



(53) 



and A*^^-* are given by, 
A(±) EE r2 T (^-3 + ri) 



(VF(1 -t-2cos6i))/(Ty) 



(55) 



It is then convenient to measure the excitation energy in 
units of T, i.e. set ti = £,pjT, for i = 1,2, 3, 4. The resulting 
(54) inelastic collision integrals reduce to 



^J.n(±) 
^ Pl 



1 



:"o(ii)V] 



pi 



[VTj 



Tin 



hl-noih)) (7r2 + i?)(^(±)(ti) / dhnoih) 



tl+t2 



sinh 



\ 2 



(±) 



(i2) 



(56) 



where I/t^^ is the quasiparticle-quasiparticle collision rate 
given in Eq. ([T]). These same transformations applied to the 
elastic collision integral in Eq. (j3Hl) yield, 

= lno(ti)(l-no(ii))vp, • (^^f) ^ X ^(±)(ii), 

(57) 

where 1 / Tei is the rate for quasiparticles on the Fermi surface 
scattering elastically off the aerogel. Finally, the left-hand 
side of the linearized transport equation provides the driving 
terms that determine the particular solution for the nonequi- 
Ibrium distribution functions. The driving terms which are 
even and odd under ti — > — ti are 

L^^Hh) = lnoih)il-no{h)) Vp, • (^^^^^) . (58) 

We simplify the even and odd components of the transport 



equation by an additional transformation of the distribution 
function. 



CW(t) 



T~ 2cosli(t/2) 



(59) 



The linearized transport equation with both inelastic and 
elastic channels included in the collision integral then reduces 
to the linear integral equationsj^^ 



1 /I 



cosh(ti/2) \ ti 



t? + 7r2 + ^)C(±)(ti 

tl - t2 



sinh((ti - t2)/2) 



C^^Ht2), (60) 



Physical solutions for ^^^^ are non- vanishing only in the low- 
energy region near the Fermi level, i.e. \t\ <^ 1. Thus, we 
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can Fourier transform, 



(61) 



and convert the integral equation for C^^' (t) into a linear 
differential equation for Q^'>{z), 



^2^2 )<?(±)(^) T 27r2A(±'sech2(^z) C,^^\z) 



-2tt sech(7rz) 



-i7rtanh(7rz) / ' 



(62) 



where 



(63) 



We cast this differential equation into standard form defined 
on the domain x G [—1,1] with the transformation, 



tanh(7rz) — x , 



and the differential operator, 



dx 



df 



dx J (1 — x^) 



(64) 
(65) 



/ • (66) 



Thus, the nonequilibrium distribution function is obtained as 
the solution of an inhomogeneous linear differential equation. 



where 



(67) 



The particular solution for ^ is obtained as an expansion, 
Z'--^ - J2 ^nMx) , (72) 



in the complete set of orthonormal eigenfunctions, {(/)„(a;)}, 
of the homogeneous differential equation. 



2?[0„] +2a„0„ = 0, 



(73) 



where a„ is the eigenvalue associated with the eigenfunction, 
(j)n{x). The coefficients, {v4„}, are determined from Eqs. 
((67)) . (|68)) and ((73)) and the orthogonality condition. 



+1 



In particular. 



where 



dx(j)*^,{x)(l)n{x) = 5„ 



+ 1 



dx<t)*^{x)R^-\x) , 



(74) 



(75) 



(76) 



is the overlap of the n"* eigenfunction with the driving term 
in the transport equation. Also, we set A^~^ = above 
and hereafter. The same term appears in the kernel of the 
heat current. Thus, the thermal conductivity, in particular, 
S'^(T), is determined by the weighted sum over the eigenvalue 
spectrum, 



Stt^ ^ - Ak 



(77) 



G. Thermal Conductivity 

The heat current for example can be expressed in terms of a 
particular solution of Eq. ((67)) , 

X rio(?p)(l-«o(Cp)) ^p(^p)- 

Carrying out the transformation from V'p(Cp) ^ Z{x) we 
obtain the following expression for the thermal conductivity, 

J(cv«/)(^'/^in)5.(T), (70) 

where 

S^{T)^^ l^\xR^-Hx)Z^-\x). (71) 



H. Pure Fermi-Liquid 

For pure '^Hc (r^j oo) Eq. ((75)) reduces to the differential 
equation for the associated Legendre functions,— 

The bounded, odd-parity solutions, relevant to heat trans- 
port, on the domain x £ [— 1, -f 1], are the associated Legen- 
dre polynomials, Pn{x) with eigenvalues, 2a„ = n{n + 1) for 
n = 1,2,3,.... The standard orthognality relation for the 
polynomials is^ 

j^\xP^{x)*PUx) = ^^^5nu'. n>l. (79) 

The evaluation of the overlap integrals leads to the solution 
for the thermal conductivity (Eq. ^) obtained by Brooker 
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and Sykes^^i^ and independently by Jensen et al^, with 
transport time = Ti„ S"^ , where 



5*: 



6 



E 



2m + 1 



TT^ ^-^^ \ [m{m + 1) — 2X^]m(m + 1) y ' 



(80) 



which depends on the angular average of the scattering am- 
plitude via A„ and is independent of temperature. Thus, the 
thermal conduetivity diverges asKocl/TasT— >0 because 
the number of thermal excitations, the heat capacity and the 
number of final states for binary collisions are all vanishing 
as T. Note that = (1^(1 + 2cos6i))/(M^) is a measure of 
the relative importance of forward vs. backscattering, and 
is restricted to the domain, — 1 < A„ < 3. The resulting 
spectral sum, S^, is finite since A^ < 3. However, at any 
fixed temperature k increases dramatically for quasipartiele 
scattering that is predominantly in the forward direction, i.e. 
for A„ 3. Note also that S°° -t^- — in this limit. 

BIT i-^A„ 



solution in the neighborhood of the singular point; in partic- 
ular, since 7 > 1 we select the physical solutions which must 
behave as 



0„ ~ {i+xr/\ x^-i 



(82) 
(83) 



Thus, we extract the behavior near the singular points and 
express 



(84) 



where gn{x) is analytic on the domain [—1, +1], and governed 
by the differential equation. 



(1 - 2:^)^ - ^^^^ + ~ 7(7+ l)].9n = . (85) 



Analytic solutions on the finite domain can be represented 
as a Taylor expansion about x = 0, 



I. Exact Solution 



Here we extend the exact solution for pure '^H o^^'^° to that 
of '^He in aerogel described by the two-channel collision in- 
tegral for binary quasipartiele collisions and quasiparticle- 
aerogel scattering. Bennett and Rice extended the analysis 
of Refs. [TTII20I to coUisional scattering of s- and d-electrons 
combined with electron-impurity scatteringJ^ Their results 
for the electrical and thermal conductivity are expressed as 
a sum over weighted integrals of products of Gegenbauer 
polynomials. Our analysis also starts from a two-channel 
extension of the integral equation of Refs. [T1II20I i.e. Eq. 
(|60)) . The results presented below are a closed form analytic 
solution to the linearized Boltzmann equation and thermal 
conductivity, and an exact perturbation theory result for the 
inelastic corrections to the elastic limit which is used to de- 
velop a very accurate approximate solution for the thermal 
conductivity that is valid for all temperatures (above Tc), 
pressures and aerogel densities, and is fast and easy to eval- 
uate. 

Elastic scattering by the aerogel modifies the form of the 
eigenfunctions for the nonequilibrium distribution function, 
and leads to an eigenvalue spectrum that varies strongly with 
temperature. The key parameter is the structure constant 
(7) in Eqs. (|62)) and (|63)) in the differential equation for 
the distribution function. The temperature dependence is 
conveniently exhibited by scaling Eq. (|63|) in terms of the 
cross-over temperature, T-,,, defined in Eq. ([5]), 




(81) 



The eigenfunctions of Eq. ([75| for any 7 > 1 must be 
bounded on the interval [—1,-1-1]. The singular points at 
a; = ±1 have indicial equations with one physically allowed 



Qnix) — ^ ^ GraX^' 



(86) 



m— 



The differential equation determines the recurrence formula 
for the coefficients. 



G 



Nr, 



m+2 



(87) 



(m + 2)(m + 1) 
iV,„ = [m{m - 1) + 2to(7 + 1) + {7(7 + 1) - 2a,„}] . 

Thus, the solutions break up into even and odd parity so- 
lutions depending on the coefficients Go and Gi. For even 
parity solutions, we set Go 7^ and generate the solutions 
from the recurrence relation: 



Gr, 



N0N2 . . . iV„_2 



Go, TO = 2,4,., 



Similarly, for the odd-parity solutions we start from Gi 7^ 
and find 



Gr, 



In either case 



NiN3...Nrr 

(m + 2)! 



Gi , m = 3, 5, . 



lim — > 1 

m-roo G,n 



(89) 



(90) 



Thus, the series solution diverges at = 1 unless the ex- 
pansion truncates at a finite value of m. This restricts the 
physical solutions for gn{x) to a set of polynomials, and an 
eigenvalue spectrum determined by the condition: 

G„+2 = ^ Ar„ = . (91) 

Expressing the eigenvalue as 2a„ = e„(e„ -t- 1), we obtain 

e„ = 7 + n, n = 0,2,4,... (n = 1,3,5,...), (92) 
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for even (odd) parity solutions. The corresponding eigen- 
functions are 



Coefficient G^""* is fixed by the normalization of (j)n{x), 



(93) 



with the summation over even (odd) integers for even (odd) 
parity eigenfunctions. The coefficients can be expressed in 
terms of Gamma functions. In particular, for the odd-parity 
eigenfunctions, which are relevant for computing the thermal 
conductivity. 



r(f+i) r(7 + f + i + ^) 
r(i^ + i) r(7 + f + i) 



(94) 



G\ 



(«) 



m—\ p— 1 



m + p + 1 



and B{x,y) = r(x)r(y)/r(a; + y) is the Beta function»i 



(95) 



We can now evaluate the spectral sum in Eq. (|77|) to obtain 
an exact solution for the thermal conductivity. In particular, 



_(7 + 2£)(7 + 2£-l)-2A« 



(96) 



where 



'2t-l 
2p-l 



.1 r(7- 



T{2p)V{i-p+l) 



(97) 



Although Eq. ((96|) provides us with an exact, closed form 
solution for the thermal conductivity over the full temper- 
ature and pressure range, < T <C T/, the sums defin- 
ing S^{T) involve ratios of Gamma functions. Thus, care 
must be taken in evaluating these functions even for moder- 
ate values of their arguments. This is particularly true in the 
low-temperature limit, T <^ T^, since the scaling parameter, 
7 becomes large. However, the limit T <^ Ti, can also be 
evaluated using perturbation theory. 



J. Perturbation Theory 



The perturbation expansion through first order becomes, 

m.,Pi) = 6i;[[dnl"'>] , (100) 
SI-[6nl"^]+Si;\[Snl'^]. (101) 







where i(^pi, Pi) represents the driving term on the left side 
of Eq. ([34|) . For heat transport the zeroth-order solution of 
Eq. (|100p is simply the distribution in the elastic scattering 
limit, 



T 



(102) 



and now provides the driving term for the first order correc- 
tion in Eq. (|10ip . This equation has the same integral kernel 
as that of Eq. pOOp and so we can express the first-order cor- 
rection in terms of an inelastic correction to the scattering 
time. 



At tempertures T <C inelastic quasiparticle collisions are 
relatively infrequent compared to clastic collisions off the 
aerogel. Thus, the inelastic collision integral in Eq. 
is of order 



6I2_ 



<5< 1, 



(98) 



and we can formally expand the integral equation and the 
deviation from local equilibrium in the small parameter S, 



(99) 



T 



(vp. • vr)Ti(ep.), 



(103) 



where ti is the first-order correction to the mean scatter- 
ing time, r^i. This distribution function gives the first order 
elastic collision integral. 



(V,, . VT) ^ 



(104) 



The solution of the Eq. (jlOip can then be expressed in terms 
of Ti. The analysis of the inelastic collision integral, eva7u- 
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ated with the zeroth order nonequihbrium distribution func- 

-(0 



tion, Sn^^ , leads to 



n =- ^(no(ii)(l-no(ti))ii)~^ 

Tin 



(105) 



where h = I{^i)/T^, with /(^i) given by Eq. [H and 



Ji = ^ y rf66no(6)i^(ei+6) 

= -^(Tr^+t^) ti(l-no(ii)) , , (106) 

with K(^) defined by Eq. 1531 The resulting first-order cor- 
rection for the collision time reduces to 

n = -i^(^^ + t?)(l-iA.) , (107) 
which vanishes in the "ballistic limit" for inelastic collisions, 

The first order correction to the thermal conductivity is cal- 
culated by evaluating Eq. (pij) with the first-order correc- 
tion, Sn^\ Writing j(^) = — (JkVT, we obtain, 



2 r + oo £2. 

5k = -Nfv] T / dt ^ Ti(t) . 

3 ^ ^ i_oo 4cosh2(t/2) ^ ^ 



(108) 




FIG. 2: Pressure dependence of the low-temperature limit of the 
thermal transport time TkT^ for pure ^He. The red are the 
data from Ref. [l^. The blue — is the spd model without en- 
forcing the FSSR, and the [black ■] are the spd model with the 
FSSR enforced. Inset: Ak calculated in the spd model with the 
FSSR enforced. 



A. Results for ^He-aerogel 



After the integration over Ti{t), and scaling to the elastic 
limit for the thermal conductivity given in Eq. p3p . we 
obtain. 




(109) 



III. RESULTS 



Theoretical models for the quasiparticle collision probabil- 
ity, W{9,(l)), for pure ^He have been proposed by a number 
of authors.— i^>2^i^ We use an extended version of the s- 
p model introduced by Dy and Pethicki^, described as the 
spd model in Sec. |Vl In Fig. [5] we compare the results for 
the thermal transport scattering time, r^T^, for pure '^He 
calculated in the spd model with the Landau parameters 
taken from Refs. Il7lfl9l and the limiting low-temperature 
thermal conductivity measurements from Greywall (Table II 
of Ref. [l6l ) . The theoretical and experimental results are in 
agreement over the pressure range, p = — 25 bar provided 
the forward-scattering sum rule (FSSR) is enforced (see Eq. 
(|165p in the Appendix). Above 25 bar there are deviations 
0% < Stk/tk, < 14% indicating the role of additional scat- 
tering not described by the spd model. Note in particular 
that the dimensionless scattering parameter, A^, is nearly 
constant over the entire pressure range, i.e. 1.0 ^ A„ < 1.3. 



Theoretical results for heat transport in '^He-aerogel based 
on the two-channel solution for the thermal conductivity are 
shown in Fig. [3l In addition to the mfp describing the aero- 
gel, the input data for bulk -^He used to generate these results 
density (n), effective mass (m*), Fermi velocity (u/) and the 
Fermi liquid parameters {F^'"') for I < 2, all of which are 
taken from the database provided in Ref. [13 and [l^. The 
Fermi-liquid parameters are used to construct the inelastic 
scattering rate using the spd model as described in the Ap- 
pendix (Sec. IV)) . 

The cross-over from the high-temperature regime dominated 
by inelastic quasiparticle collisions to the low-temperature 
regime dominated by elastic scattering by the disordered 
medium occurs over a fairly broad temperature range for 
dilute aerogels with long mfp. The elastic regime below 
T ~ 5 - lOmK is well described by k = k,, + 0{T^), with 
K^i given by Eq. ([55]) . The pressure dependence of the slope 
of Kei(T), while not visible in Fig. [31 is shown clearly in Fig. 
m Note that limT->o can provide a determination of the 
clastic mfp for the aerogel. 

Figurc[4]for k/T highlights the deviations in the thermal con- 
ductivity from the elastic limit limit even at temperatures of 
order a few milli-Kelvin. Similarly, in the high temperature 
limit the product, kT, approaches the bulk '^He limit de- 
termined by inelastic scattering. Significant deviations from 
the pure '^He limit arc shown in Fig. Oover a wide range of 
temperatures above T*. 
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FIG. 3: Theoretical results for the thermal conductivity of He- 
aerogel vs. T and p for an elastic mfp of £ = 1700 A are shown in 
black — • — ; results for pure ''He are shown as green lines. The 
pressure ranges from p = — 32 bar in steps of 2 bar starting with 
the the upper curve at p = bar. 



k/T [erg/mK^ -cm-sec] 




ln(T[mK]) 

FIG. 4: k/T vs. ln(r[mK]) and p[bar]. The inputs and labels are 
the as Fig. [3] The slopes for the elastic limit, k^i/T, are shown 
in darkred . 



B. Comparison with Experiments 



Barker et.aS- reported results for the thermal conductivity of 
^He in 98% at aerogel of k = 7.2 mW/mK at p = 32.4 bar 
and T = 2.20 mK. They also added two monolayers for 
*He which displaces the solid ^He coating the silica aero- 
gel strands, and measured a slight increase in the thermal 



FIG. 5: kT vs. r[mK] and p[bar]. The inputs and labels are the 
same as those of Fig. |3] Note that results for kT for bulk "'He are 
shown as darkgreen — . 



conductivity, i.e. k = 7.7 mW/mK at T = 2.22 mK. Com- 
parison of these two data points with the theoretical predic- 
tions for this pressure arc shown in Fig. [6] 

If the difference in the mfp with and without the ^He is at- 
tributable to spin-exchange scattering of itinerant '^He spins 
by the localized solid ^He spins,— then we can estimate 
the contribution to the scattering rate from indirect spin- 
exchange scattering to be, 



1 



Vf 



1 



1 



(110) 



and thus a mean time for spin-exchange scattering of T,pi„ ~ 
0.15 /.tscc, i.e. several orders of magnitude longer than the 
mean time for elastic scattering off the aerogel strands, r^i ~ 
€3ho+4ho/^^/ — 8.6 ns. For scattering off a random distribution 
of Ns localized spins via a Kondo interaction, 

u = - ^( J.„,/n)S, • cr 5(r - R,;) , (111) 

i=i 

the Born approximation implies an additional contribution 
to the scattering rate. 



{J^^JnY SiS + 1) . 



(112) 



Thus. we estimate the indirect exchange interaction to be 



4 

9^ 



Ef 



1/2 



0.5 mK/spin , 



(113) 

which is in agreement with the order of magnitude estimate 
for Ji„d inferred from the absence of a low-field Ai — A2 tran- 
sition in ■^Hc-aerogel<^ 
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■-1.0 -0.5 0,0 0.5 1.0 1.5 2.0 2.5 3.0 



lii(T [mK]) 

FIG. 6: Comparison k/T vs. In(T[mK]) with experiments. Data 
shown as ■ is from the Lancaster group for a 98% aerogel and 
a pressure of p ~ Obar.— Data shown as ■ (purple ■) is from 
the Stanford group with two monolayers of ^He coating (without 
''He) the silica strands at high pressure (p — 32.4 bar), also for 
98% aerogel but grown in a different laboratory. The — show 
limr^o k/T. 



The Lancaster group also reports results for the low- 
temperature (i.e. T <^ T^^) thermal conductivity of normal 
■^He-aerogel at low pressures, for aerogels with porosities of 
95% and 98%3 Results for g = 98% reported in Ref. 
yield a much smaller mfp, £ = 950 A, than the Stanford data, 
suggesting significant differences in aerogels of the same den- 
sity prepared under different growth conditions. Note that 
the authors of Ref. [3^ attribute the deviations from the the- 
oretical curve onsetting near T ~ 20 mK (ln(T[mK]) ~ 3) 
to Kaptiza boundary conductance through the experimental 
cell walls. 

Although these results provide estimates for the aerogel mfp 
they do not provide a test of the the theory. Measurements of 
the thermal conductivity over the full temperature and pres- 
sure range of normal ■^He-aerogel should provide a strong test 
of the two-channel theory based on homogeneous disorder 
since we have an exact solution for the thermal conductiv- 
ity in this model. Conversely, if significant deviations from 
the theoretical predictions are observed they could indicate 
new physics associated scattering and transport of fermionic 
excitations in a correlated random medium. 



C. Scaling Function 

The exact solution for the thermal conductivity in the two- 
channel scattering theory for '^He-aerogel can be expressed 
in terms of a scaling function. Normalizing k by the thermal 
conductivity in the elastic limit, k^i, from Eq. p3p gives, 



— = ^x5.(T) = F(T/T*,A0 



(114) 



Note that t^i/t^^ = {T/T*Y, and that S^{T) calculated from 
Eq. ((77)) [Eq. ([M|) in Sec. Ill I| provides the exact scaling 
function, F{x, A), since S^{T) depends only on x = T/T* and 
the scattering ratio, —1 < < 3. Thus, the test of the two- 
channel transport theory would be to demonstrate that the 
thermal conductivity of ^ He-aerogel obeys the scaling behavior 
over the full temperature and pressure range of the normal 
state, and a broad range of aerogel density and mfp. 

The exact solution for k/k^i x {T/Ti,) = xF{x,Xi^) is shown 
in Fig. (O. The calculation of spectral sum, S^{T), was car- 
ried out using arbitrary-precision floating point arithmetic in 
order to evaluate the ratios of the Gamma functions or large 
arguments that enter Eq. (j96[) with sufficient precision to 
obtain accurate results for the triple sum that defines (T) . 
In particular, the points labeled "exact" in Fig. ([7]) were ob- 
tained with the floating point precision set at 55 digits and 



each sum was cutoff after 30 terms were computed. One can 
obtain reasonably good results with a lower precision setting 
for the floating point arithmetic, but double precision on a 
32-bit machine limits the accuracy of the results, particularly 
in the limit T < T^,. 

Also shown in Fig. ([7]) are calculations of the scaling function 
based on an approximate analytic formula that is numerically 
fast and easy to evaluate. The approximate scaling function 
is constructed from the the asymptotic limits for S^{T) for 
T ':^Ti, and T 0, as well as the leading order perturbative 
result for T <^Ti,, as described below. 

The limiting behavior for the exact scaling function is known 
from the asymptotic limit, a; ^ 1, and perturbation theory 
about X = 0. In particular. 



F{x,K) = l 1 



where 5*" is given by Eq. (jSO)) . 



l-f7r2(l-iA«)x-2 ,a;«l 



,a; > 1 



(115) 



The most common approximate solution for multi-channel 
scattering is based Matthiessen's Rule, which in this context 
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T/T* 

FIG. 7: Scaling function for the thermal conductivity plotted as xF{x,\k) = k/k^i x (T/Ti,) as a function of a; = T/T-^ for Ak = 
0.0, 0.5, 1.0, 1.5, 2.0, 2.5, starting from the lowest to the highest curve, respectively. The exact results based on Eqs. I81l97l and l96l are 
shown as the Black — • — , the Cyan ■ are based on the approximate analytic formula given in Eqs. (| 1321 11331 1 1181 and 1 127p . and the 
Orange - are the results for pure, bulk ^He normalized to the elastic limit for "'He-aerogel. 
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FIG. 8: Comparison of several approximates with the exact scal- 
ing function, x F{x,Xk), plotted as a function of a; = T/T* for 
Are = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, starting from the lowest to the 
highest curve, respectively. The exact results are shown in Black 
— • — , while the approximate scaling functions, labelled as poly, 
exp and MR, are described in the text. 



can be expressed as 



1 

rMR 



(116) 



i.e. the total transport scattering rate is the sum of indepen- 
dent rates for purely elastic and purely bulk inelastic trans- 
port. The resulting expression for the the thermal conduc- 
tivity, normalized to its value in the elastic limit, 



^ = FMH(^/^^A,), 



(117) 



defines the approximate scaling function, i^„n(x, A„) given 
by, 



1 



l-(S'~)-ia;2 ,2;<1 



1 + (S^)-^ 



,a; > 1 . 
(118) 



The MR scaling function deviates from the exact result of 
Eq. (|115|) for the leading order finite temperature correction. 
Curiously, the exact result for the leading order correction 
is equal to that obtained from Fmb.{x,Xk,) by approximating 
S'^ with just the first term of the sum in Eq. ([50]) . This 
approximation is very good in the limit of nearly forward 
scattering. In this limit the inelastic channel leads to large 
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thermal transport for T >Ti,. As a result the MR scaling 
function gives a very good approximation to the exact scal- 
ing function in the limit of large for all x. This is shown 
clearly in Fig. ([5|). However, the MR scaling function devi- 
ates from the exact scaling function when backscattering in 
the inelastic channel is significant, i.e. for A„ < 1.0. These 
deviations are also clearly visible in Fig. ([8]). 

We can try to improve on the MR scaling function by incor- 
porating the exact perturbative result for F{x^ A^) for x <C 1. 
We construct an interpolation formula that connects the ex- 
act asymptotic limits. A simple extension of Matthiessen's 
interpolation formula is the two-parameter, rational polyno- 
mial function, 



1 



1 



which has the limiting forms, 



1 - (a - l)x2 ,a;<l 

- + - — , a; > 1 . 

a b x'^ 



(119) 



(120) 



We then fix the coefhcients from the exact asymptotic limits 
for F(x,Ak) in Eq. (fTTS]) . 



aS"- - 1 



(121) 
(122) 



Although this approximate scaling function works well for 
the .T <C 1; it docs a poor job in the intermediate and high- 
temperature region a; > 1 (green curves in Fig. |8]), and is 
particularly poor for A^ — ^ 3. If we consider the leading 
order correction to the asymptotic limit a; — > oo we obtain 



coo 1 1 

i^^^+C74-+0(-). 

x^ x^ x^ 

For the polynomial approximate we obtain, 

^4 ' = "(1+ 6|i(l_Aj3))2 ' 

while the MR scaling function gives 

cr = -(^r)'- 



(123) 



(124) 



(125) 



Both approximate scaling functions give the correct sign for 
the leading order correction, however in the limit A^ — > 3, 
where we know the MR scaling function approaches the exact 
result, we see that Cf"' is large and negative. 



4 



5 
6^ 



1 



(1-AJ3)2 



(126) 



whereas C^^^^ 



— 1. This discrepancy in i^poiy is traced to 
the contamination of the temperature region a; > 1 by the 
exact solution that is valid for a; ^ 1. 



We might remedy this problem with a two-parameter inter- 
polation that limits the contamination between a; ^ 1 and 
a; 1. In particular, consider the approximate scaling func- 
tion. 



-e 
2 



For a; < 1, 



> 1 - aa;^ + 0(a;'') . 



(127) 



(128) 



Note that there arc only exponentially small corrections to 
the leading order result for a; ^ 1 coming from the terms that 
are fixed by the asymptotic solution for a; 3> 1. Similarly, for 
a; ^ 1, the term that is fixed by the exact solution for a: <C 1 
is now exponentially small and we obtain, 

Using these expansions and the exact leading order asymp- 
totic limits we obtain 



a 

b = S?° 



^(1-AJ3) 



(130) 
(131) 

This two-parameter interpolation formula yields a better ap- 
proximation to the exact scaling function, particularly for 
Ak < 1. However, F^^^ under estimates the maximum in 
X F{x, Xk), and this deviation is enhanced as A^ — >■ 3, as is 
clear from Fig. ([5]). The basic result of this analysis is that 
the MR scaling function, Fmr , is accurate in the limit of large 
Ak, but deviates from exact scaling for A^ ^ 1- By contrast 
the two-parameter exponential scaling function, F^.^^, is ac- 
curate in limit A^ < 1, but shows increasing errors from exact 
scaling in the cross-over region, a; ^ 0{1), for 1 < A„ < 3. 
This suggests that we combine these two scaling functions 
into a single scaling function by weighting the respective re- 
gions of accurate scaling, i.e. 

^lpprox(a:, A«) = p(A^)_F„p(a:, A„) 

+ (l-p(A«)) FMR(a;,AJ, (132) 

where the weight function p(Ak;) is chosen on the physical 
domain, —1 < A^ < 3, to satisfy, p(— 1) = 1, p(+3) = 0. 
Thus, the simplest weight functions which map the physical 
domain onto the interval [0, 1] are 



P(A«) = 



1 + A« 



(133) 



The quadratic weight function, i.e. s = 2, leads to remark- 
ably good agreement with the exact scaling function for the 
entire domain of A^ and reduced temperature, x ~ T/T^,. 
This comparison is shown in Fig. ((T]). Note that the maxi- 
mum deviation for any of the computed values is less than 
0.6%, and careful examination shows that these small errors 
occur near the maxima of x F(x, A^)- Thus, the main result 
here is that Eqs. 11181 11271 11321 and 11331 provide numerically 
fast and accurate formulas for calculating the thermal con- 
ductivity over the full temperature and pressure range within 
the two-channel scattering theory for normal "^He-aerogel. 



C Scaling Function 
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D. Scaling for ^He-aerogel 

The analysis of the pressure dependence of the thermal con- 
ductivity of pure ^He based on the spd scattering amplitude 
described in Sec. IIIII and App. |V] implies that the thermal 
transport scattering parameter is nearly pressure indepen- 
dent, i.e. Ak — 1.3 for 5 bar < p < 34 bar with a smooth 
drop to Ak — 1-0 as pressures between 5 and Obar (see inset 
of Fig. El). 

Pressure independence of the scattering parameter, A^, im- 
plies that the thermal conductivity for all temperatures above 
the superfluid transition, all pressures and all elastic mean- 
free paths should collapse to a single scaling function when 
normalized to its value in the elastic scattering limit, i.e. 
liuiT^o K = 'foi given in Eq. Thus, for ^He-aerogel we 

expect that thermal conductivity for all T, p and £ to col- 
lapse to the narrow band of scaling functions shown in Fig. 
ini A complete set of measurements of the thermal conduc- 
tivity of •^He-aerogel for all T, p and a wide range of aerogel 
mfp would provide a strong test of this theory, particularly 
the assumption of uncorrelated disorder described by a single 
mfp. 



inates, to a low temperature regime dominated by elastic 
scattering off the aerogel. The exact solution to the two- 
channel Boltzmann-Landau transport equation reported here 
provides quantitative predictions for heat transport in liq- 
uid '^He-aerogel. An approximate solution derived from the 
asymptotic solutions and perturbation theory is accurate to 
less than 0.6%. A key result of this work is the scaling 
function, F(T/T*,A„), that describes the exact solution for 
the normalized thermal conductivity, k/k^i, for all pressures, 
temperatures (above Tc) and aerogel density. A complete set 
of measurements of the thermal conductivity of ■^He-aerogel 
for all T, p and a wide range of aerogel mfp would pro- 
vide a strong test of this theory, particularly the predicted 
scaling behavior based on two-channel scattering and the 
assumption of homogeneous disorder described by a single 
mfp. Conversely, systematic deviations from the predicted 
scaling function behavior should provide a quantitative mea- 
sure of the role of fractal correlations associated with the 
structure of the aerogel. The limited data that is available 
already hints that two-channel scattering is insufficient and 
that spin-exchange scattering between itinerant '^He spins 
and localized '^Hc spins contributes to the low-temperature 
thermal conductivity. 
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FIG. 9: Scaling of the thermal conductivity with T/T* for ^He- 
aerogel for the physical range: 1.0 < < 1.3. Inset: Scaling 
function for a much larger parameter range of transport scatter- 
ing parameters, including Ak = 0.0, 2.0, 2.5, in addition to the 
physical range. 



IV. SUMMARY 



Liquid '^He impregnated into silica aerogel is a model sys- 
tem for investigating the effects of quenched disorder on 
the properties of a strongly correlated Fermi liquid. In the 
normal Fermi liquid the transport of heat, mass and spin 
by fcrmionic excitations exhibits cross-over behavior from 
a high temperature regime, where inelastic scattering dom- 
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APPENDIX: ^HE SCATTERING AMPLITUDE 



The binary collision amplitude for quasiparticles in pure '^He 
in the low-energy region near the Fermi surface depends on 
the momenta and the spin state of the initial and final pair 
of excitations. In particular the dimensionless scattering am- 
plitude is 



P3"3 

Taia2;a3a4(PliP2;P3,P4) = 2iV/ X 

PlCkl 



P4Q;4 




P2a2 (134) 



where t is formally defined by the matrix elements of a tran- 
sition operator between incoming (1,2) and outgoing (3,4) 
quasiparticles. For a Fermi liquid with only exchange inter- 
actions such as ■^He the total spin S and any one component, 
Sz, are conserved by collisions. For an unpolarized Fermi 
liquid there is no preferred direction for the spins to align. 
As a result all three spin-triplet amplitudes are equal and 
there are only two independent amplitudes corresponding to 
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the total spin S = and S = 1, which we label as the singlet 
(s) and triplet (t) amplitudes, 



Ts = - [Tti;n - Tf4.4t - T^lf.n + Tit;it] (135) 
1 , , 



^tt;n = Ui;U ■ 



(136) 



Also note that amplitudes which differ by inversion of all the 
spin projections are equal, 

Thus, we use a short-hand notation, 







-i—l'i ■J--J- 


(138) 




= Tn;n 




(139) 








(140) 



and express the spin-projection amplitudes in terms of the 
singlet and triplet amplitudes 



iff — it 

Tn = J(Tt + T,) 



(141) 
(142) 

(143) 



The T-matrix can then be expressed in terms of Tj and the 
corresponding symmetric (triplet) and anti-symmetric (sin- 
glet) spin matrix elements, 



(-) 



(144) 



where 



^Qia2;a3Q4 ~ g ('^"l"3'^a2a4 ^ '^QiQ4 '^a2a3 ) ■ (1^5) 

Since there are only two independent amplitudes it is often 
useful to use the symmetric and anti- symmetric amplitudes 
defined as 

= ^(Ttt + Tti) = ^(STt + T,) (146) 



- ^ (Tn - Tti) = J (Tt - T,) 



Inverting, we have 



Tt = + T°) 
Ts = (T"' - 31") . 



(147) 



(148) 
(149) 



The two sets of amplitudes, T*'° or Tt,s, define different, but 
equivalent representations for the spin-dependent T matrix. 
The T'''" amplitudes are the amplitudes for the T-matrix ex- 
pressed in terms of the direct "particle- hole" channel, 1 — > 3 
and 2-^4, 



^ QlCt2;CK3Ct4 



(150) 



For quasiparticle scattering on the Fermi surface the scatter- 
ing amplitudes, Tt,s, reduce to functions of the directions of 
the quasiparticle momenta on the Fermi surface. 



Tt,s(Pl,P2;P3,P4) Tt,s(pi,P2;P3,P4) 



(151) 



Furthermore, rotational invariance implies that Tt^^ can be 
expressed in terms of the relative direction cosines. 



X2 = P2 • Pl = cos U = P3 • P4 , 
X3 = P3 • Pl = COS 6*3 = P4 • P2 , 
X4 = P4 ■ Pl = COS 6*4 = P3 • P3 . 



(152) 
(153) 
(154) 



The fourth column of equalities follows from momentum con- 
servation for I Pi I = Pf, 

Pl+P2=P3 + P4- (155) 

The conservation law also implies that there are only two 
independent angles. We adopt Abrikosov and Khalatnikov's 
parametrizationi in terms of the angle 6 between the two 
incoming momenta, and ^, the angle between the planes de- 
fined by n = Pl X p2 and n' = p3 x p4. 



n • n 



COS( 



|n||n'| 

Thus, Tt,,(pi,p2;p3,P4) =Tt^,(6l, 



3^3 - Xj 
1 — X2 



(156) 



The Pauli exclusion principle requires the T-matrix to be 
anti-symmetric under exchange of either the initial or the fi- 
nal state of the two fermions. Thus, the spin-singlet (triplet) 
amplitude is necessarily symmetric (anti-symmetric) under 
exchange of the initial or final momenta, or in terms of the 
scattering angle. 



(157) 
(158) 



Thus, we can formally expand the singlet (triplet) amplitudes 
as a sum over even (odd) functions of cos(to0), 

T,(6',(/)) = 4™)(cos6') cos(m0), (159) 



ni—0 



(160) 



m—l 



Note that Tt vanishes for cj) — Tr/2 and = 3tt/2. For these 
angles the momentum transfer in the direct and exchange 
channels is identical, in which case exchange symmetry re- 
quires the triplet amplitude to vanish identically. 

Microscopic analysis of the two-particle propagator and its 
relation to the quasiparticle scattering amplitude leads to 
an identity between the scattering amplitude in the forward 
direction, and the Landau parameters, F^'°', that define the 
quasiparticle molecular fields. In terms of the symmetric 
and anti-symmetric amplitudes in the p-h channel. Landau's 
identity for the forward scattering amplitude is^^. 



,4> = 0) 



E 

t=0 



(161) 
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where 



F 



l + F/'7(2£+l) ■ 
In terms of the singlet and triplet amplitudes for 
Ts{0,<p = 0) = ^{AI~3A''e)Ve{cos0 

Tt(^^,0 = O) = ^(A,^+A?)7'Hcos^) 



(162) 

= 0, 
(163) 

(164) 



Exchange symmetry leads to an additional constraint on the 
triplet amplitude. In the limit 6 ~ 0, (f> — 0, i.e. for pi = 
P2 = P3 = P4i the triplet amplitude necessarily vanishes. 
Thus, from Eq. (|164p we obtain the forward scattering sum 
rule (FSSR), 



lim TtiO,0) = yiA', 



Af) = 0. 



(165) 



s-p-d Scattering 



Several microscopic and phenomenological theories have 
been proposed for the quasiparticle scattering amplitude in 
3iiq13^i24^i26^ Here we adopt a slightly modified version of 
the model proposed by Dy and PethickJ^ They proposed 
a minimal model for the scattering amplitude that obeys 
exchange anti-symmetry. In particular, if we assume the 
singlet and triplet scattering amplitudes are to a good ap- 
proximation given by the m = and m = 1 terms, we have 
Ts ~ As (cos 6*) and Tt ~ At (cos 6') cosi^. In this case we can 
fix the expansion coefficients of At^s{cos6) in terms of the 
forward-scattering amplitudes, A^'", and thus the Landau 
parameters, Fp'""^ 



Tt ~ ^ (A| + A'l) Viicos 9) cos ( 



(166) 
(167) 



l>0 



The quasiparticle lifetime, Ti„(T) in Eq. [U as well as the 
thermal transport time, Tk,{T), due to binary quasiparticle 
collisions in pure '^He are determined by angular averages of 
the spin-averaged transition probability. 



1 

Tin 



N' 



f 



VfPf 

5r(A.)ri„ 



{W) [ksTf 



with 



where W = 
in Eq. EH 



A« = {W{l + 2cose))/{W), 



(168) 
(169) 

(170) 



W-fi and the angular average is defined 



probability can be expressed in terms of the dimensionlcss 
singlet and triplet scattering amplitudes. 



W = iT^r + 3|Ttr -l-25R{TtT:} 
and the quasiparticle lifetime becomes. 



(171) 



(172) 



Note that for weighted averages of W in which the weight 
function is even under exchange (i.e. (f> ^ + tt) the cross 
term in Eq. (jl7ip vanishes. Similarly, for the thermal trans- 
port time we can write = A^/ (W) where 



A. 



((H-2cos6») W). 



(173) 



In the spd model the Fermi-surface average of the rate be- 
comes. 



(W) = {AD + -{A',) 



(174) 



We evaluate this rate in terms of Legendre expansion of the 
forward-scattering amplitudes. For either singlet or triplet 
channel, 



A(cos6') =yAePi cos W) 

i 

The angular average of A^ is given by 

{AD = ' 



where 

= [ dxPi{2x^ - l)Pe'{2x^ - 1) . 
Jo 

Similarly, for the angular averages of the form, 
{A^ (1 + 2cos0)) ^ ' 

with 



(175) 



(176) 



(177) 



(178) 



Writing Wab = ^h^^N jWab, the transition 



Lu'= [ dx {Ax^ -l)Pe{2x^ -l)Pt{2x^ -1). (179) 
Jo 

These coefficients are listed in Table U for £,£' < 2. 

The input for our calculations of the transport properties of 
bulk '^He as well as •^Hc-aerogel are the Fermi-liquid parame- 
ters. The measured values of these parameters are collected 
in Ref. [l^, and are also available onlineJ^ The Landau in- 
teraction parameters, Fq, Ff , and Fq are accurately known 
from measurements of the heat capacity, first-sound veloc- 
ity and magnetic susceptibility of pure normal '^He, while 
determinations of F| and F^ have also been obtained from 
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the spd model, or the accuracy of the known Fermi liquid 
data is inadequate, or both. The most problematic aspect of 
the spd model as it stands is that the FSSR is badly violated, 
when evaluated with y4|'° = for £ > 3. In particular the 
largest violation in the FSSR, 



TABLE I: Coefficients defining the angular averages of (VK) (Cui) 
and {W{1 + 2cosS)) (Lui) in the spd scattering model. 



Serror 



(180) 



measurements of the zero sound velocity and spin-wave res- 
onance for normal "^He, respectively. However, these param- 
eters are not as accurately determined. Less is known about 
the magnitude and pressure dependence of the £ = 2 contri- 
bution to the exchange interaction, .F^r^ and much less is 
known quantitatively about the Landau interaction parame- 
ters corresponding to harmonics £ > 2, although evidence of 
interactions in higher order scattering channels is suggested 
by the observation of high frequency pair exciton^ modes in 
superfluid ^He-B.— 

The scattering model we use throughout is defined by Eqs. 
(|166p and (|167p with the added assumption that we truncate 
the expansion, i.e. set = for ^ > 3. This approxima- 
tion is reasonable if the contributions to the Fermi-surface 
averages of the scattering rate fall off sufficiently rapidly with 
increasing £ > 2. 

The spd model with the Fermi liquid data for £ < 2 as input 
qualitatively describes the decrease in the transport time, 
TkT'^, with increasing pressure (shown in Fig. [2]) and is 
within 25% of the experimental values for r^T^ over the full 
pressure range. However, the comparison clearly shows that 



is at low pressures, which is also where the discrepancy (refer 
to Fig. [2]) between theory (solid line) and experiment (red 
diamonds) is greatest. This is a significant violation of the 
Pauli exclusion principle, and is an indication that either the 
determinations of ^-rc inaccurate, that there is significant 
weight in the interaction channels for £ > 3, or both. 

Respecting the Pauli exclusion principle, by enforcing the 
FSSR, is likely more important than knowing precisely the 
distribution of higher angular momentum channels that ac- 
count for the missing weight in Eq. 11801 Thus, we enforce 
the FSSR by fixing the least known material parameter in 
the spd model, i.e. we replace 



^2 -^^2 



^2 



(181) 



The importance of enforcing the FSSR appears to be born 
out by the improvement between theory (black dots) and ex- 
periment shown in Fig. ^ For pressures below p < 25 bar the 
agreement is nearly perfect. Thus, the deviations between 
theory and experiment at higher pressures likely reflects real 
limitations of the spd model, i.e. there is scattering that re- 
duces heat transport that is outside the spd scattering model. 
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